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A study of the in-plane resistance noise on a high quality single crystal of Lai.97Sro.o3Cu04 deep 
inside the spin-glass phase reveals the onset of non-Gaussianity and the insensitivity of the noise 
to the in-plane magnetic field. The results indicate that the charge dynamics becomes increasingly 
slow and correlated as temperature T ^ 0. The analysis of the higher order noise statistics provides 
evidence for the existence of a collective ground state of charge clusters, which seem to coexist with 
charge-poor antiferromagnetic domains that are frozen at such low T. 
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I. INTRODUCTION 

There is growing experimental evidence for the emer- 
gence of spatial inhomogeneities at the nanoscale in 
many strongly correlated electron systems. In addition 
to charge, other degrees of freedom, such as spin and 
lattice, typically also play an important role in these 
materials, leading to the competition of several ground 
states and the resulting nanoscale phase separation. ^"'^ 
Such nanophases have been variously described as bub- 
bles, stripes, clumps, clusters, or domains, depending 
on system details. Other investigations have discussed 
the emergence of inhomogeneous charge-ordered phases 
in analogy with the smectic and nematic states in liq- 
uid crystals.^ In general, many different configurations of 
such nanoscopic ordered regions often have comparable 
(free) energies, such that these metastable states are sep- 
arated by barriers with a wide distribution of heights and, 
thus, relaxation times. This leads to the slow dynamics 
typical of glassy or out-of-equilibrium systems. In fact, 
the frustration caused by the competition of interactions 
on different length scales may give rise to the emergence 
of an exponentially large number of metastable configura- 
tions and the associated glassy dynamics even in the ab- 
sence of disorder This scenario is relevant to many ma- 
terials, including cuprates, where a long-range Coulomb 
repulsion competes with a short-range attraction that re- 
sults from magnetic exchange interactions. The disorder 
should further stabilize the glassy phase. Even though 
the emergence of glassiness thus appears to be almost 
ubiquitous at low temperatures (T), the glassy charge 
dynamics and out-of-equilibrium behavior in general re- 
main poorly understood. 

In cuprates, various experimental techniques provide 
evidence for glassiness in the spin sector at T < Tsg{x), 
where Tsg is the spin glass transition temperature and 
X is the doping)^"— In La2Sr2_a;Cu04 (LSCO), where 
the long-range antiferromagnetic (AF) order of the par- 



ent compound La2Cu04 is completely suppressed for 
X w 0.02 but two-dimensional short-range AF correla- 
tions persist^, the high-temperature superconductivity 
(HTSC) sets in at a; « 0.055. The spin glass (SG) phase 
emerges with the first added holes and extends all the 
way to slightly overdoped x ~ 0.19 (Refs. [rllTof). thus 
coexisting with HTSC. Unlike conventional spin glasses, 
here the low-T phase results from cooperative freezing of 
moments in different domains in which spins are antifer- 
romagnetically orderedJ^ Therefore, the spin glass has a 
local "stripe" character and is often called a "cluster spin 
glass". On the other hand, other experiments suggest 
that charge is clustered in antiphase boundarie o^'^'^' '— 
that separate the hole-poor AF domains in Cu02 (ab) 
planesi^ii^"— Imaging techniques on other cuprates (see, 
e.g., Ref. [20l ) have confirmed the presence of randomly 
distributed charge domains with short-range order in the 
underdoped or pseudogap regime. It is plausible to ex- 
pect, based on general arguments, that the ground state 
of such a random distribution of nanoscale charge do- 
mains is a charge glass. The nature of the ground state is 
clearly of great interest, because it is from this state that 
the HTSC emerges with hole doping. However, studies 
of charge glassiness, i.e. of dynamic as opposed to static 
charge ordering, have been relatively scarce.^1 Moreover, 
experiments at relatively high T may be difficult to inter- 
pret in many materials because of the changes in struc- 
tural or magnetic symmetry. At low T, on the other 
hand, the charge dynamics may become so slow that the 
charge distribution appears to be static on experimental 
time scales. 

Therefore, the resistance noise spectroscopy employed 
at very low T and on very long time (t) scales is clearly 
a technique that is well suited for probing the charge 
dynamics. It has already been used successfully in 
charge (Coulomb) glasses in nonmagnetic systems i^^Ti^i 
The analysis of the noise statistics can provide evidence 
for the presence of a large number of metastable states, 
which are an essential prerequisite for well-known dy- 
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FIG. 1: (Color online) The temperature dependence of the 
zero- field cooled in-plane resistance Rab below 30 K (Ref. [53). 
The solid line is a fit with slope To- Inset: Schematic phase 
diagram of hole-doped LSCO. 



namical features of glassiness, such as aging, memory, 
and breaking of ergodicity. 

In LSCO, the noise has been measured in the c-axis 
or out-of-plane resistance Rc of the x = 0.03 sample)^ 
which is located in the pseudogap regime, but it does 
not superconduct at any T (Fig. [T] inset). The data 
have provided evidence that, deep within the SG phase 
at T <iC TsG where transport is insulating, the charge 
dynamics become increasingly slow and correlated, i.e. 
glassy, as T — i> (Ref. Ilsl) . Moreover, the results strongly 
suggest that doped holes form a cluster charge glass. 
These conclusions are further supported by impedance 
spectroscopy.— 

In the same regime where Rc noise shows evidence 
for charge glassiness, both c-axis^^^ and in-plane 
magnetoresistance-^ (MR) exhibit the emergence of a 
strong, positive component for all orientations of the 
magnetic field B. This positive MR (pMR) shows a hys- 
teresis and memory, indicative of the presence of corre- 
lated magnetic domains. A careful and comprehensive 
analysis of the datapS^i^ points to the picture of AF do- 
mains that are frozen at low T ^ Tsg, and holes con- 
fined to the domain "walls" {i.e. patchy areas separat- 
ing the AF domains). While low B produces some mo- 
tion of the domain walls leading to a hysteretic MR, the 
main transport mechanism that gives rise to the pMR 
remains unchanged and, in fact, is the same as that ob- 
served in various nonmagnetic, disordered materials with 
strong Coulomb interactions. Basically, in a system that 
conducts via variable-range hopping (VRH) , the Zeeman 
splitting in the presence of a Coulomb repulsion between 
two holes in the same disorder-localized state leads to 
a pMR by blocking certain hopping channelsi ^^'^^ Much 
higher B leads to the reorientation of the weak ferro- 
magnetic (FM) moments2Si2i associated with each AF 
domain and oriented along the c axis, resulting in a nega- 
tive MR.^ Therefore, since the magnetic background for 
a fixed B is frozen at T <g; Tsg, at least on experimental 
time scales, only holes in the domain walls contribute to 
transport and glassiness in the resistance noise. For this 



reason, the statistics of the Rc noise, which reflects a col- 
lective wandering of the system among many metastable 
states, is unaffected by the presence of a fixed B, and it 
is also independent of the magnetic history. 

While the previous study of the Rc noise in 
Lai.97Sro.o3Cu04 strongly suggests^^ that the doped 
holes form a cluster glass state as a result of Coulomb 
interactions, albeit in the presence of a frozen, random 
magnetic background, the noise in the in-plane transport 
has not been investigated. In addition, only the effect of 
-B II c was explored. However, the in-plane transport 
is generally agreed to be more relevant to the physics of 
cuprates and, in particular, to the emergence of HTSC at 
higher charge-carrier concentrations. Furthermore, vari- 
ous theoretical models, such as those that consider fluc- 
tuations of the local electronic order and the effect on 
transport noise^^^ have focused on the in-plane trans- 
port and in-plane symmetry breaking. It should be noted 
that the relationship between the in-plane and out-of- 
plane transport in these highly anisotropic materials is 
far from trivial, and their sometimes contrasting behav- 
ior has been one of the most unusual characteristics of 
the cuprates. Indeed, comparative studies of intralayer 
and interlayer properties continue to provide key insights 
into the physics of these strongly correlated systems 

Therefore, here we investigate the noise in the in-plane 
resistance Rat of Lai.97Sro.o3Cu04 at very low T ^ Tsg 
and with B \\ ah. Just like what was observed^ for the 
c-axis transport and a different field orientation {B \\ c), 
the fluctuations of the in-planc resistance Rab with time 
provide evidence for the increasingly slow, correlated dy- 
namics as T — >■ 0. Here, however, the magnitude of the 
resistance noise is significantly smaller and the emergence 
of charge glassiness takes place at even lower T than for 
the out-of-plane transport. In spite of these quantita- 
tive differences, the overall data strongly support earlier 
conclusions^^"^'' and provide evidence that doped holes 
in lightly doped LSCO form a collective, glassy state of 
charge domains or clusters located in Cu02 planes. 



II. EXPERIMENT 

A high quality single crystal of LSCO with x — 
0.03 was grown by the traveling-solvent floating-zone 
technique. Detailed measurements were performed on 
a sample that was cut out along the main crystallo- 
graphic axes and polished into a bar with dimensions 
2.10 X 0.44 X 0.42 mm^, suitable for direct Rab mea- 
surements. Electrical contacts were made by evaporating 
gold on polished crystal surfaces, followed by annealing 
in air at 700°C. For current contacts, the whole area of 
two opposing side faces was covered with gold to ensure 
a uniform current flow through the sample. In turn, the 
voltage contacts were made narrow (~ 80 ^m) in order 
to minimize the uncertainty in the absolute values of the 
resistance. Gold leads were attached to the sample using 
Dupont 6838 silver paste. This was followed by the heat 
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FIG. 2: (Color online) Schematic of the circuit for the five- 
point ac bridge measurement of the voltage (resistance) fluc- 
tuations. Ballast resistors Ri and R2 (~ 100 MQ) balance 
the bridge, thus minimizing the effects of temperature and 
current fluctuations. Two PAR124A lock-in amplifiers detect 
the in-phase and out-of-phase components of the signal. The 
power spectrum of the former contains both sample resistance 
noise and background, while the power spectrum of the latter 
contains only the background noise. The total current was 
measured with an ITHACO 1211 current preamplifier and a 
PAR124A lock-in amplifier. 



treatment at 450°C in the flow of oxygen for 6 minutes. 
The resulting contact resistances were less than 1 £7 at 
room T. 

The sample resistance was measured using the stan- 
dard four-probe ac technique (typically at ~ 7 Hz) in 
the Ohmic regime, at T down to 0.05 K realized in a 
dilution refrigerator. The low- frequency noise was mea- 
sured using a conventional five-probe ac bridge method^ 
with the bridge current / « 100 nA in the Ohmic regime 
(Fig. [5]). This method makes the resistance noise mea- 
surement nearly insensitive to both T and current fluc- 
tuations. For example, it was confirmed experimentally 
that the correlation between T and voltage fluctuations 
decreases and becomes negligible as the balance of the 
bridge is increased. Likewise, the correlation between 
the fluctuations of voltage and current was also negligi- 
ble. Two lock-in amplifiers 7 Hz) were used to detect 
the difference voltage. This measurement setup makes 
it possible to remove the contribution of the background 
noise, which is due to the Johnson noise and noise from 
the preamplifieri^ In this experiment, the power spec- 
trum of the background noise was always negligibly small 
and white. (The background noise measured with zero 
current was white for all T and B and 2-3 orders of mag- 
nitude smaller than the sample resistance noise.) The 
output filters of the lock-in amplifiers served as an an- 
tialiasing device. Most of the noise spectra were obtained 
in the / = 10~* — 10~^ Hz bandwidth, where the upper 
bound was set by the low frequency of /, limited by the 
low cutoff frequency of RC filters used to reduce external 
electromagnetic noise as well as by the resistance of the 
sample. 

The sample is located in the SG region of the phase di- 



agram (Fig. [U inset) atT <TsG^ 7 -8K {Rei.M). In 
that T regime, Rab shows insulating behavior. In particu- 
lar, it obeys the VRH dependence Rab = Ro s^p{To/T)'^, 
To ~ 108 K, with an exponent /i = 1/3 all the way up to 
T = 30 K (Fig. [J). {Rab is henceforth denoted by R for 
simplicity.) The value /x = 1/3 is characteristic of two- 
dimensional (2D) VRH and it is in agreement with early 
results on ceramic LSCO samples where it was shown 
that this exponent is doping dependent. The localiza- 
tion length obtained from the VRH fit^ is ^ ~ 90 A, 
which is somewhat larger than the average magnetic cor- 
relation length 40 A at this doping.Siii The most prob- 
able hopping distanceii rh{T) ~ ^{Tq/TY thus spans on 
the order of ten AF domains at low T. Therefore, for the 
VRH process, the system appears to be uniform. 



III. RESISTANCE NOISE 
A. Noise in zero magnetic field 

Figure (Sl^a) shows the typical time traces of the rela- 
tive changes in resistance AR{t)/{R{t)) in zero field at 
various T. It is obvious that, at the lowest T [see also 
Fig. Elb)], the system exhibits random fiuctuations on 
many different time scales, including slow changes over 
several hours. While the Rc noise occasionally exhibited 
also some discrete switching events;^ such fluctuations 
have not been observed in the in-plane resistance noise 
measurements. On the other hand, similar to the c-axis 
transport, here the magnitude of the noise does not ap- 
pear to depend on T, but the character of the noise is 
significantly altered with increasing T. 

In particular, at the lowest T [e.g. Fig. [3l^b)], the 
noise is not Gaussian. This is clearly seen already in 
the histograms of AR{t) / {R{t)) values or the probability 
density function (PDF) of the fiuctuations. In general, 
the low-T PDF has a complex structure with more than 
one peak and, moreover, its precise shape depends ran- 
domly on the observation time. This is demonstrated, 
for example, by the PDFs of the noise data obtained 
during four sequential 3-hour intervals at the same low T 
[Fig. Hl^a)]. The wandering of the PDF shape with time 
indicates that different states contribute to the resistance 
as a function of time, signaling that the system is non- 
ergodic (glassy) on experimental time scales. Although 
the PDF broadens when the sampling time is increased to 
12 hours, since the system has more time to explore the 
free energy landscape, the PDF remains non-Gaussian 
[Fig. lUb)]. When T is increased, on the other hand, 
non-Gaussian effects become less pronounced [Fig.HJJb)], 
and vanish completely already at T ~ 0.2 K where PDF 
is described by the Gaussian distribution [Figure |4l[c)]. 

The normalized power spectra 5'i?(/) (where / is the 
frequency) of the noise AR{t) / (R) obey the well-known 
empirical law Sr oc 1//" [Fig. O^a)]. The exponent a 
increases as T is reduced [Fig. Ob)] , indicating a shift of 
the spectral weight toward lower / and the slowing down 
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FIG. 3: (Color online) (a) Relative fluctuations of the in-plane 
resistance AR{t) / (R) = (R— {R))/{R) as a function of time 
for B = at several T. (R) is the time-averaged resistance 
determined from a four-probe measurement at each T. All 
traces are shifted vertically for clarity, (b) The T — 0.064 K 
data in (a) shown on an expanded scale for clarity. 
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FIG. 4: (Color online) (a) Probability density function (PDF) 
vs. AR{t)/{R) for four sequential 3- hour time intervals at 
T = 0.082 K and in B = 0. (b) PDF vs. AR{t)/{R) for 
12-hour intervals at several T with B — 0. (c) PDF values 
at r = 0.183 K fitted to a Gaussian distribution, as shown 
{S = 5.6 X 10"^). 



B. Noise in magnetic field 



of the dynamics. A careful analysis of the noise mag- 
nitude Sii{f) at fixed / also shows that slow dynamic 
contributions {i.e. low /) to conductivity become in- 
creasingly important with decreasing T [Fig.[5jc)]. More- 
over, large values of the exponent a at low T (a > 1, up 
to « 1.5) reflect the increasing non-Gaussianity of the 
noise as T — ^ 0. Similar behavior was observed in the Re 
noisej25. as well as in some other charge22i2^ and spin'*^''*'^ 
glasses. However, we note that, while the in-plane noise 
becomes nearly white {i.e. a w 0) already at T ~ 0.2 K, 
the out-of-plane noise retains its 1// frequency depen- 
dence {a 1) even at T ^ 0.3 Ki^ This difference, 
together with the absence (presence) of switching events 
in the in-plane (out-of-plane) noise, suggests, not too sur- 
prisingly, the existence of some additional processes that 
contribute to the c-axis transport and noise. Indeed, the 
in-plane value of the power spectral density is approxi- 
mately two orders of magnitude smaller than the out-of- 
plane value 



The relative changes of R with time have been mea- 
sured also in magnetic fields B || a6 of up to 9 T at a 
fixed, low T [Fig.[6Ua)], where the zero- field noise is non- 
Gaussian. Similar to the c-axis case^ the data do not 
seem to show any effect of the magnetic field on either 
the amplitude or the character of the noise. For exam- 
ple, the low-T PDF remains non-Gaussian [Fig. ^h)] up 
to the maximum applied field. The noise power spectra 
Sr oc 1//" do not show any effect of B either [Fig.[7]Ja)]: 
the large low-T values of the exponent a do not vary with 
B [Fig. [Tjb)], and the magnitude of the noise is indepen- 
dent of B in the entire experimental range of frequencies 
[Fig. El^c)]. Furthermore, while the resistance itself de- 
pends strongly on the cooling procedure^^ all the noise 
characteristics are independent of the magnetic history. 
In particular, the same results are obtained after field- 
cooling and zero-field cooling protocols, again in analogy 
with Ref. m. 
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FIG. 5: (Color online) (a) Octave-averaged power spectra 
S'_r(/) shown for different T and -6 = have been corrected 
for the white background noise. Solid lines are linear least- 
squares fits to the form Sr oc 1//". (b) The exponent a vs. 
T. The dashed line guides the eye. (c) S'h(/) for / = 0.1, 1.0, 
and 7.2 mHz, determined from the fits in (a), as a function of 
T. The error bars are standard deviations of the data. As T is 
reduced, (high-) low-frequency contributions to conductivity 
become (decreasingly) increasingly important. 



C. Higher-order statistics 

Many noise sources in nature are Gaussian. For Gaus- 
sian random processes in general, all n-point (n > 2) 
correlation functions can be expressed in terms of prod- 
ucts of pairwise correlation functions.—!^ Therefore, the 
measurement of higher-order correlations cannot give any 
new information beyond that which is already contained 
in the pairwise correlation function (the Fourier trans- 
form of the power spectrum). 

A superposition of many independent sources gives rise 
to a Gaussian noise. For example, noninteracting, iden- 
tical sources or fluctuators (e.g. electron emission from 
a cathode) produce noise that is random, not correlated 
in time, so that the corresponding power spectrum is in- 
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FIG. 6: (Color onhne) (a) Noise AR{t)/{R) vs. time at T = 
0.118 K for several B \\ ah. All traces are shifted vertically 
for clarity, (b) PDF vs. AR/R in different fields for 12-hour 
observation periods; T = 0.118 K. 



dependent of frequency ( "white" ) . On the other hand, a 
superposition of many noninteracting two-state systems 
(TSSs) with a broad distribution of characteristic times 
gives rise to a Gaussian noise with a power spectrum of 
1/f type. In such a material, the spectrum S at each 
/ is given by the sum of fixed components from differ- 
ent TSS. Repeated measurements should thus yield the 
same values of the power spectrum at a given /, within 
experimental accuracy. In reality, since actual measure- 
ments take place over finite times, repeatedly measured 
S fluctuates randomly around its mean value, which is 
time- independent. Thus the power spectrum of the time 
series of such measurements of S (the "second spectrum" ) 
is white; that is, the noise is Gaussian. 

In many complex systems, however, strong deviations 
of the noise statistics from Gaussianity have been ob- 
served, indicating that the noise cannot be described by 
noninteracting models4^^— For example, repeated mea- 
surements reveal the "wandering" of the power spectrum 
at a given / with time, rather than settling into a sum 
of fixed components coming from different fluctuators 4^ 
This wandering means that, during the time scale of each 
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FIG. 7: (Color online) (a) The octave-averaged power spectra 
S'_r(/), shown for several B at T = 0.118 K, have been cor- 
rected for the white background noise. Solid lines are linear 
least-squares fits to the form Sr oc 1//" . (b) The exponent a 
vs. B at T = 0.118 K. (c) Si{(/) for / = 0.1, 1.0 and 7.2 mHz, 
determined from the fits in (a), vs. B ai T — 0.118 K. The 
error bars are standard deviations of the data. 



measurement, the system "visits" only a small number of 
states in the free energy landscape {e.g. several shallow 
valleys within one deep valley). In that case, the value 
of 5 at a given / is determined by the characteristic re- 
laxation times of the shallow valleys. After a sufhciently 
long time, the system undergoes a transition to another 
deep valley with a different set of relaxation times for 
transitions between its shallow valleys. Repeated mea- 
surements of 5 at a given / will thus yield values that 
are correlated in time, with the corresponding nonwhite 
second spectrum. This type of system is obviously non- 
ergodic. 

The use of higher-order correlation and response func- 
tions is now broadly accepted as a standard tool for ana- 
lyzing dynamical heterogeneities in a wide variety of sys- 
tems. In particular, the second spectrum, a fourth-order 
noise statistic, was first used to investigate the dynamics 
of spin glasses and provide information about correla- 
tions among local fluctuators JSsH The second spectrum 



S2if2, /) is defined'*^ as the power spectrum of the fluctu- 
ations of Siiif) with time and, in practice, is calculated 
in narrow frequency bands {e.g. octaves) / — {fL,fH)- 
It was introduced as a technique for understanding how 
noise power redistributes among narrow frequency bands 
as a system evolves in time. The second spectrum's util- 
ity is based on its sensitivity to correlations among the 
fluctuators that give rise to the non-Gaussian proper- 
ties of the noise. As discussed above, if the noise arises 
from an ensemble of statistically independent fluctuators 
(Gaussian noise), S'2(/2,/) is white {i.e. independent of 
/2). On the other hand, S'2(/2,/) « 17/2 '^ for interact- 
ing fluctuators.— i^i^ 

The second spectra have been calculated using digital 
filteringi^i^ First, the fast Fourier transform (FFT) is 
applied to the data, which are then digitally bandpass 
flltered in a given frequency range / = (/l,///) with 
/h — 2/l. This allows the analysis of AR{t)/ (i?) in any 
bandpass. After taking the inverse FFT of the flltered 
data, they are squared point by point to obtain the vari- 
ance of the noise vs. time for a given /. Finally, the power 
spectrum of the time-dependent variance of the noise sig- 
nal {i.e. the second spectrum) is calculated. 52(/2,/) is 
then divided by the square of the averaged total power 
over unit time in the bandwidth limited SR{f) to obtain 
the normalized second spectrum. Figure IHa) illustrates 
the results of that analysis. The normalized second spec- 
tra, with the Gaussian background subtracted, are shown 
for the lowest and highest T in the experiment. The ex- 
ponent (1 — /3), which is a measure of correlations, clearly 
increases with decreasing temperature. Therefore, as T is 
reduced, the fluctuators become strongly correlated and 
the non-Gaussian nature of the noise becomes more pro- 
nounced. A detailed dependence of (1 — /3) on both T 
and B is given in Fig. [SJb) . Just like other noise char- 
acteristics, (1 — P) exhibits i?-independent behavior at a 
flxed T, indicating that the magnetic field does not affect 
the nature of the correlations. 

The second spectrum has been useful also in distin- 
guishing between different kinetic models that have iden- 
tical spectra. In particular, the scaling of ^2 with respect 
to / and /2 has been used in studies of spin^i^ and 
Coulomb glasses^^Ti^i to discriminate between two rival 
theoretical approaches - the interacting droplet model 
and a picture of hierarchical dynamics. In the droplet 
approachf^Si^ which relies crucially on the finite range 
of the interactions and assumes compact droplets, the 
slow events {i.e. low-/ noise) originate from rearrange- 
ments ("flipping") of a smaller number of large droplets, 
while the higher-/ events come from a larger number of 
smaller elements that are faster. "'^■'^^ Furthermore, ac- 
cording to this picture, large droplets are more likely to 
interact than small ones with a probability depending on 
their actual size, leading to the conclusion that S2 will be 
stronger for lower /. However, second spectra obtained 
for different / need to be compared for a flxed /2// {i.e. 
on time scales determined by the timescales of the fluctu- 
ations being measured), since spectra taken over a flxed 
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FIG. 8: (Color online) (a) Normalized second spectra S'2(/2), 
with the Gaussian background subtracted, for two T in B = 0. 
Solid lines are fits to S2 oc 1/ f2~^ ■ (b) The exponent (f — /3) 
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data. The dashed line guides the eye. The results shown in 
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and 2-4 mHz octaves. 
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FIG. 9: (Color online) Second spectrum S'2(/2,/) measured 
in octaves / = (/l,2/l) (a) in zero magnetic field at T = 
0.082 K and (b) in B = 6 T at T = 0.118 K. Dashed lines are 
linear fits to guide the eye. No dependence on the magnetic 
field has been observed within the scatter of data. 



IV. CONCLUSIONS 



time interval average the high-frequency data more than 
the low-frequency data. Hence, in the interacting droplet 
model, 5*2 (/2, /) should be a decreasing function of / at 
constant /2//. In the hierarchical picture of glasses, 
on the other hand, there is no characteristic scale: the 
second spectra depend only on the ratio /2//, not on /, 
and thus can be collapsed onto a single curve. Such scale 
invariant 6*2 (/2,/), consistent with the hierarchical pic- 
ture of glassy dynamics, have been reported in systems 
with long-range interactions, such as conventional spin 
glasses {e.g. CuMn)i^i^ and 2D Coulomb glassesi^"— 

Our findings for the second spectra obtained for differ- 
ent values of / in both zero and finite magnetic fields are 
shown in Figs.[ni[a) and[Hl[b), respectively. It is apparent 
that ^2 decreases with increasing / for a fixed /2//, con- 
sistent with generalized models of interacting, compact 
droplets or clusters. Therefore, this result further sup- 
ports the picture of spatial segregation of holes into inter- 
acting, hole-rich droplets or clusters, which are separated 
by hole-poor AF domains. These dynamic charge inho- 
mogeneities thus result from the competition of short- 
range interactions, revealed by S'2(/2,/) (Fig-©, and a 
long-range Coulomb interaction. 



Our measurements of the in-plane resistance noise in 
X = 0.03 LSCO have demonstrated the emergence of 
slow, correlated dynamics and nonergodic behavior at 
very low T <^ Tsg- The gradual enhancement of this 
glassy behavior with decreasing temperature strongly 
suggests that the phase transition to a charge glass state 
occurs at T = 0. Several findings in our study can be used 
to rule out spins as the origin of the observed glasslike 
behavior. For example, the resistance noise spectroscopy 
reveals that the significant arrest of the charge motion 
in Cu02 planes occurs only at T < 0.2 K {i.e. deep 
inside the spin-glass phase). Furthermore, while both 
in-plane and out-of-plane resistance are strongly affected 
by the magnetic field^^i^l all the noise characteristics are 
insensitive to both the magnetic field and the magnetic 
history, further indicating that the observed glassiness 
reflects the dynamics of charge, not spins. This is analo- 
gous to the magnetic insensitivity of the noise in a glassy, 
fully spin-polarized 2D electron systemi^ We note that in 
conventional spin glasses, by contrast, all the resistance 
noise characteristics are affected by j3 i42,43,53 

Therefore, this study provides evidence for the ex- 
istence of a cluster charge glass ground state in a; = 
0.03 LSCO, in agreement with the conclusions based 
on noise measurements on c-axis oriented samples^, 
dielectric^^ , and magnetotransport^ studies. The strong 
non-Gaussianity of the noise indicates that the dynam- 



8 



ics of charge clusters cannot be attributed merely to the 
switching between frozen and non-frozen states of indi- 
vidual clusters (TSS) with different characteristics, but 
rather that the charge clusters are correlated. Although 
they affect the interlayer or c-axis transport )^^'^^ here we 
have shown that the correlated charge clusters are lo- 
cated in Cu02 planes, where they seem to coexist with 
hole-poor AF domains that remain frozen at very low 
temperatures T ^ Tsg, at least on experimental time 
scales. Further work is needed to examine the evolution 
of this dynamically inhomogeneous state with doping and 



its possible coexistence with HTSC. 
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